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by Marcus F. Heldmann and Paul R. Wieber 
Lewis Research Center 

SUMMARY 

Vaporization histories of n-heptane drops in a combustor with superimposed travel- 
ing transverse oscillations were computed by using drop- evaporation theory. Oscilla- 
tions in the rates of vaporization of an array of repetitively injected drops in the combus- 
tor were obtained from summations of individual drop histories. A frequency response 
of the entire vaporization process to oscillations in pressure was evaluated. The typical 
response curve showed a maximum value at a specific frequency. The response factor 
approached zero at lower frequencies and a constant negative value at higher frequencies, 
which indicated a potential for both the driving and damping of acoustic oscillations. 

Changes in combustor and drop parameters were investigated. The critical fre- 
quency, which divided regions of positive and negative response, increased with an in- 
crease in combustor pressure, combustor axial gas velocity, and amplitude of the pres- 
sure oscillations. It decreased with an increase in drop radius. Changes in initial drop 
velocity and temperature caused no significant change in critical frequency. An evalua- 
tion of liquid-oxygen vaporization showed propellant properties to be important. 

A normalized response characteristic based on a correlation of vaporization parame- 
ters is presented. The correlation shows that the peak response occurs when a drop is 
vaporized in a time period approximately equal to the oscillation period. Negative gains 
occur when the vaporization time is greater than about three times the oscillation period. 
Drop lifetimes less than three times the oscillation period are frequently encountered in 
unstable rocket combustors. 


Part of the material given herein was presented at the Second Combustion Insta- 
bility Conference, Los Angeles, California, Nov. 1-5, 1965. 


INTRODUCTION 


System analyses of unstable combustion (refs. 1 and 2) have shown that propellant 
vaporization is a process which, when varied, can control system stability. In such sys- 
tem analyses, the dynamic response of the rate of propellant vaporization to the acoustic 
pressure oscillations is not readily identified. It is the purpose of this report to evaluate 
the response characteristic for a range of conditions in vaporization parameters. The 
evaluations provide a measure of the sensitivity of the vaporization process to acoustic 
oscillations by indicating its driving or damping potential in an unstable feedback loop. 

The drop vaporization theory, developed in reference 3 to evaluate combustor perfor- 
mance for a vaporization- lirnited combustion process, is used in this study to compute 
drop vaporization histories in a rocket combustor with superimposed acoustic oscillations. 
In the vaporization calculations, the acoustic oscillations affect the drop acceleration 
and heat and mass transfer processes by giving the drop three-dimensional velocity com- 
ponents and causing perturbations in drop temperature and evaporation rate. With these 
variations, the rate of vaporization depends on the frequency of the oscillation, and a fre- 
quency response can be evaluated. This application of the theory differs from that used 
in the system analyses of references 1 and 2, where the vaporization rate varied only 
with perturbations in drop Reynolds number about a mean condition in vaporization and 
was independent of the frequency of the oscillations. 

Calculations are made for n-heptane drops vaporizing in a cylindrical combustor 
containing heptane-oxygen combustion products. Pressure, velocity, and temperature 
oscillations, associated with the first traveling transverse acoustic mode, are super- 
imposed on the normal combustion flow process. Individual drop histories are obtained 
for boundary conditions which include changes in combustor pressure, final combustion 
gas velocity, drop radius, initial drop velocity, initial drop temperature, drop radial 
position in the combustor, and the amplitude and frequency of the acoustic mode. An 
evaluation of oxygen vaporization is made to explore the effects of propellant properties 
on vaporization in an acoustic field. 

Individual drop histories depend on the time and position of drop injection into the 
combustor with respect to the orientation of the acoustic mode. A summation of these 
individual histories is used to obtain the time variations in the vaporization rate of re- 
petitively injected drops throughout the combustor. The summation provides a response 
of the entire vaporization process treated as an entity. This response factor is ex- 
pressed as the ratio of the percentage oscillation in vaporization rate to the percentage 
oscillation in pressure and accounts for phase relations. Variations in the response fac- 
tor with frequency and boundary condition are presented. 
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SYMBOLS 


d 

a 

C P 

D v 

F 

f 

H v 

J i 

k 

L 

& 

m f 


m 

N 

n 

P c 

P L 

Pr 

P T 

AP c 

R 

R d 

R d, 0 
Re 

% 


contraction ratio, dimensionless 
speed of sound, ft/sec 
specific heat, Btu/(lb)(°R) 
diffusivity, ft^/sec 


transformed frequency, 



, cps 


frequency, cps 

heat of vaporization, Btu/lb 

Bessel function of order i 

thermal conductivity, Btu/(ft)(sec)(°R) 

length of combustion zone, ft 

burning rate parameter, dimensionless 

final Mach number, Up/a, dimensionless 

molecular weight, lb/lb mole 

fraction of weight burned or vaporized per unit length, ft” ^ 
response factor, (AW/W)/(AP c /P c ), dimensionless 
interaction index (see ref. 6) 

2 

average gas pressure, lb/in. 

2 

vapor pressure at liquid surface, lb/ft 
Prandtl number, (c p ^ m M m ) A m 

o 

total gas pressure, 144(P c + P ), lb/n 

2 

maximum peak-to-peak pressure oscillation in combustor, lb/in. 

gas constant, (ft)(lb)/(lb mole)(°R) 

drop radius, ft 

initial drop radius, microns 

Reynolds number, (2R ( j P m AV )/M m 

combustor radius, ft 

drop radial position, ft 
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Sc Schmidt number, |U /D n 

’ nr v^m 

T c gas temperature associated with average gas pressure, °R 
drop temperature, °R 
T T total gas temperature, T c + T c , °R 
t time, sec 

t^ cycle time, 1/f, sec 
U axial gas velocity, ft/sec 
U d axial drop velocity, ft/sec ' 

Up final axial gas velocity, ft/sec 

V transverse gas velocity, ft/sec 

V d transverse drop velocity, ft/ sec 

AV vector sum of all directional components of velocity difference between gas and 
drop, [(U - U d ) 2 + (V r - V d> r ) 2 + (V Q - V d) 0 ) 2 ] 1/2 , ft/sec 
W vaporization rate of drop array, lb/sec 
w drop weight, lb 

w vaporization rate of single drop, lb/sec 

y ratio of specific heats, dimensionless 

9 drop angular position, radians 

ju viscosity, lb/(ft) (sec) 

p density, lb/ft^ 

t characteristic time (see ref. 6), sec 

Subscripts: 

L liquid value 

m mean value based on film composition and/or temperature 
peak frequency at which response factor reaches maximum 
r radial component 

v vapor value 

6 angular component 

0 initial value 
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50 value when 50 percent of drop mass is evaporated 
Superscripts: 

~ perturbation component 

rate with respect to time 
— average 


METHOD OF ANALYSIS 

The technique used herein is described in two parts. The first part presents the as- 
sumptions in the model and lists the equations which are solved to give single drop histo- 
ries. By means of a sample calculation, the second part depicts the use of single drop 
histories in generating the vaporization response factor for a column of repetitively in- 
jected drops exposed to an acoustic oscillation. 


Vaporization Model 

Drops of n-heptane are assumed to be vaporizing in combustion gases, composed of 
stoichiometric reaction products with oxygen in a cylindrical combustor, with an estab- 
lished traveling-transverse resonant acoustic mode. The pressure and gas velocity os- 
cillations attributed to the acoustic mode are derived from the relations in references 4 
and 5: 
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The wave is assumed to be adiabatic with 
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( 5 ) 


The oscillations are superimposed on the mean level of the parameters affecting drop 
evaporation and motion. In the theory of reference 3 used in this study, vaporization rate 
is controlled by heat and mass transfer at the drop surface. Drop motion is obtained 
from a momentum balance in three vector directions. Axial gas velocity is proportional 
to the fraction of drop mass vaporized, and it attains a final assumed velocity at complete 
vaporization. 

The drop history is described by following equations for weight evaporation rate, 
heating rate, and acceleration in an axial direction. The weight evaporation rate can be 
defined as 
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(6) 


The heating rate T L is found by 


T r = 
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2 + 0. 6 Sc 1/3 Re 1//2 
2+0.6 Pr 1 / 3 Re 1//2 


( 8 ) 


The final equation necessary for describing the drop history is that of acceleration in an 
axial direction, which is 
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( 9 ) 


Similar equations can be written for the Cartesian components of the two transverse 
drop accelerations V d r and V d q, where the quantity U - U d is replaced by the trans- 
verse relative velocity between the gas and the drop and transformed into corresponding 
Cartesian coordinates. The axial gas velocity is given by 


u = u F 



( 10 ) 


Liquid and gas properties were taken from equations in references 3 and 6 plus curve 
fits of data for heat of vaporization, vapor pressure, and liquid density. Direct analytical 
solutions are not possible, and evaluations are obtained by numerical techniques which 
employ a high-speed digital computer. 


Vaporization Histories 

The dynamic behavior of the entire vaporization process is obtained from a summa- 
tion of instantaneous values in the vaporization histories of individual drops. This is best 
described by the procedure for a specific set of boundary conditions. The example 
chosen is for drops of 50-micron radius injected at a velocity of 100 feet per second and 
an initial temperature of 650° R. Mean combustor pressure is 300 pounds per square 
inch, and the final gas velocity attains a value of 800 feet per second. The acoustic os- 
cillation has a peak-to-peak pressure amplitude of 20 percent of combustor pressure at a 
frequency of 3000 cps. Drop vaporization histories with these boundary conditions depend 
on two variables: (1) the phase time of drop injection with respect to the oscillation in 
pressure, and (2) the radial position of injection within the cross section of the combustor. 
Since the acoustic mode has traveling wave properties, these two variables are sufficient 
to specify the vaporization of any drop injected into the combustor. 

The vaporization histories of an array of continuously injected drops at a radial posi- 
tion r/R^ of 0. 912 will be considered first. Figure 1 shows the vaporization history 
of a drop injected at the mean pressure preceding a rise in oscillatory pressure 
(27rft -0 = 0 in eqs. (1) to (3)). Vaporization rates, in addition to other combustor pa- 
rameters, are shown as a function of cycle time. For comparison, the vaporization his- 
tory without acoustic oscillations is also shown. Although the oscillation substantially 
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Figure 1. - Time histories of drop and combustor parame- 
ters in steady-state and oscillating flow field. Drop 
conditions: radius, 50 microns; velocity, 100 feet per 
second; temperature, 650° R; radial position, 0.912 of 
combustor radius. Gas conditions: pressure, 300 pounds 
per square inch; final velocity, 800 feet per second; oscil- 
lation frequency, 3000 cps ; oscillation peak-to-peak 
amplitude, 0.2 of average gas pressure. 


reduces the total time for vaporization, the characteristic of interest in this study is the 
oscillation in vaporization rate. Oscillations in vaporization rate occur at twice the fre- 
quency of the pressure oscillations. This rectified-type response is caused by the sensi- 
tivity of the vaporization process to the absolute velocity difference between the drop and 
combustion gases. Velocity difference reaches a maximum twice during each pressure 
oscillation and never becomes equal to zero because the two transverse gas velocity com- 
ponents have a 90° phase difference. 


Summation of Histories 

A specific summation of single-drop histories is used to evaluate the time variations 
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Drop vaporization rate, 


l 


Summation times 



Cycle time, rad 

Figure 2. - Generation of array of drops from individual 
drop histories. Drop conditions: radius, 50 microns; 
velocity, 100 feet per second; temperature, 650° R; 
radial position, 0.912 of combustor radius. Gas condi- 
tions: pressure, 300 pounds per square inch; final 
velocity, 800 feet per second; oscillation frequency, 
3000 cps; oscillation peak-to-peak amplitude, 0.2 of 
average gas pressure. Each vaporization curve begins 
at zero point. 


i 


in vaporization rate of a one- dimensional ar- 
ray of repetitively injected drops. The drop 
histories in figure 2 illustrate the summa- 
tion. 

Figure 2 shows vaporization rate histo- 
ries for drops injected at a fixed radial and 
angular position in the combustor. Drop his- 
tories displayed along the ordinate are re- 
lated to the cycle time of the acoustic oscilla- 
tion along the abscissa. The pressure oscil- 
lation is shown along the abscissa for refer- 
ence. 

Four drops are injected every cycle at 
the four cycle times of 0, jt/2, n, and 37 t/ 2. 
Vaporization histories vary among drops in- 
jected at different times during one pressure 
oscillation, as shown by the first four histo- 
ries. However, drops injected at times 2n 
apart experience identical acoustic pressure 
and velocity fields and, thus, have identical 
histories. Eventually, the same number of 
drops are completely vaporized per cycle as 
are injected, and the number of drops in the 
array becomes constant over each full cycle. 
For the example illustrated, the fully devel- 
oped array contains 15 drops, decreasing in 
mass down the chamber and ranging in age 
from a new drop just injected to an old one 
almost completely vaporized. With the array 
fully developed, the distribution of vaporiza- 
tion rates along the array can be obtained at 
any time during the oscillation in pressure. 
This is done in the example at cycle times 
of 7r , 3 tt/ 2, 0, and 7 t/ 2, as indicated by the 
vertical dashed lines of figure 2. 

The mass vaporizing from the entire ar- 
ray of drops at each specific cycle time is 
the summation of that vaporizing from all 
the individual drops. 
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W = 2 Wi (11) 

i 


This summation is plotted as a function of cycle 
time in figure 3. The curve of figure 3 is the re- 
sponse of the vaporization rate of the array to an 
acoustic oscillation. The average of all the sum- 
mations over a complete cycle is also shown. 

The effect of angular and radial displacement of 
the drop in the chamber due to the oscillating flow 
field was small and was neglected in relating the 
phase of vaporization rate and pressure. The 
vaporization rate is higher throughout the array 
at both the maximum and minimum pressure con- 
dition in the oscillation than at the mean pressure 
conditions. Lower rates occur at mean pressure because the total velocity is a minimum 
at the mean pressure in a traveling wave. 

The example was simplified for clarity. In actual computation, the response curve 
was improved in two ways. First, the minimum number of drop histories that would give 
an acceptable summation was set at 16 to improve the average value. For the boundary 
conditions of the example case, this was obtained by injecting eight drops per cycle in- 
stead of four. Thus, the number of drops injected per cycle was dependent on the drop 
lifetime and the period of the oscillation. Second, the continuity of the curve was im- 
proved by sampling the array at 16 cycle times for all cases instead of only 4 times as 
shown. 
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Cycle time, rad 

Figure 3. - Vaporization rates of array of con- 
tinuously injected drops during oscillation in 
pressure. Drop conditions: radius, 

50 microns; velocity, 100 feet per second; 
temperature, 650° R ; radial position, 0.912 
combustor radius. Gas conditions: pres- 
sure, 300 pounds per square inch; final 
velocity, 800 feet per second; oscillation 
frequency, 3000 cps; oscillation peak-to- 
peak amplitude, 0.2 of average gas pressure. 


Response Factor 

The response curve of figure 3 shows that a sinusoidal oscillation in pressure causes 
an oscillation in vaporization rate which exhibits harmonic components of the basic fre- 
quency. The relation between pressure and vaporization rate, therefore, cannot be ex- 
pressed simply by a gain and phase angle, as is frequently done in dynamic analysis. 

A response factor (effective gain), which combines the amplitude ratio and phase re- 
lation of vaporization rate and oscillation pressure in a single quantity, was numerically 
evaluated. This response factor is assumed positive when vaporization rate and chamber 
pressure are both above or below their mean values and assumed negative when vaporiza- 
tion rate and chamber pressure are on opposite sides of their means. This is applied 
throughout the cycle. The regions where the response factors are positive and negative 
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are indicated in figure 3. When this convention is used, the response factor for the en- 
tire oscillation can be expressed by the ratio 

AW 

N = -55L (12) 

AP c 



The vaporization term is expressed by the summation 

£(w-w) + 

AW _ J) 

W 277 

0 

and the pressure term is defined by a similar summation, evaluated in all cases for the 
maximum pressure amplitudes that occur at the combustor wall. 

The response factor is an index of the degree of driving or damping which the vapor- 
ization process provides in a feedback loop, as generally postulated for combustion insta- 
bility (ref. 7). It is comparable to the pressure exponent in a linearization of 

W~P» 

where 

dW 


dP c 


The response factor in this example is -0. 177. Vaporization rate contributes to 
negative feedback when rate and pressure are interrelated. The value applies to all an- 
gular positions at 0. 912 of the maximum radius. It is independent of angular position be- 
cause of the traveling wave characteristics of this transverse mode. 


(14) 


(15) 
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TABLE I. - RANGE OF BOUNDARY CONDITIONS 


[Combustor temperature, T ' 5000° R. ] 


Final axial 
gas velocity, 
U 



Initial axial 
drop velocity, 

U d, 0’ 
ft/sec 

a 100 


50 

200 

a 100 

4oo 


Initial drop 
temperature, 

T L, 0’ 


*650 


1 

500 

800 


Boundary conditions of reference case. 


Area Average 

The response factor varies with ra- 
dial position because of the changes in 
pressure and velocity oscillations speci- 
fied by equations (1) to (3). Pressure os- 
cillations for this mode of resonance van- 
ish at the axis of the combustor, although 
velocity oscillations persist throughout 
the combustor. A response factor, which 
accounts for this effect, is obtained from 
an average of the response factors at three 
radial positions. The three radial posi- 
tions are the mean radii of three concen- 
tric equal areas. These radii are 0. 912, 

0. 707, and 0. 409 of the maximum radius. 


Range of Boundary Conditions 

By using the ideas just presented, the response of the entire vaporization process is 
evaluated for the range of boundary conditions shown in table I. These boundary condi- 
tions survey the effect of combustor pressure, final gas velocity, initial drop radius, ve- 
locity, temperature, radial injection position, and the amplitude of the pressure oscilla- 
tion. Pressure amplitude ratio is given at the circumference of the combustor, where it 
is a maximum. 

For table I, calculations were made for all combinations (approximately 600) of the 
following: 

f 

Drop radius, R d Q , microns ’50,150,500 

Pressure amplitude ratio, AP C /P C 0.1, '0. 2, 0. 4 

Radial position, r /R^, 0- 409, 0. 707, 0. 912 

Frequency, f, cps 100, 300, 500, 1000, 3000, 10 000 

The effect of propellant properties was explored by evaluating several cases of 
liquid-oxygen vaporization at selected boundary conditions. 


t 


Boundary conditions of reference case. 


12 



RESULTS AND DISCUSSION 


The variation in response factor with frequency is defined as the frequency response 
of the vaporization process. The effect of boundary conditions, illustrated in figures 4 
to 10, was obtained from a parametric study of the reference boundary conditions, as in- 
dicated in table I. Reference boundary conditions include the following: drop radius, 

50 microns; combustor pressure, 300 pounds per square inch; pressure amplitude ratio, 
0. 2; final axial gas velocity, 800 feet per second; initial axial drop velocity, 100 feet per 
second; and initial drop temperature, 650° R. 

Radial Position 

The frequency response for the reference boundary conditions is shown in figure 4. 
The response at the three radial positions and the average for the entire vaporization 
process are shown. These frequency response characteristics are typical for all condi- 
tions. The response factor exhibits a positive maximum, or peak, value at a specific 
frequency. It approaches zero at a lower frequency and decreases to a negative value at 
higher frequencies. This behavior is characteristic of a damped resonant system, where 
a peak response occurs at some frequency and a negative feedback occurs at a higher fre- 
quency. If this behavior is realistic, it is significant to the problem of combustion insta- 
bility in rocket combustors. It implies that the vaporization process may be tuned to the 
acoustic frequency of the combustion cavity. 

In general, larger positive and negative re- 
sponse factors were observed as radial position 
approached the combustor wall. The frequency 
for peak or maximum response, however, is rela- 
tively unaffected. The average response factor 
for the entire vaporization process is approximate- 
ly the same as that for the 0. 707 radial position, 
which divides the combustor into equal concentric 
areas. 

A maximum response factor of 0. 76 is obtain- 
able at a radial position of 0. 912. In a resonant 
system, where acoustic energy may accumulate, 
self- sustained oscillations are possible with a re- 
sponse factor greater than zero (the exact value 
depending on losses). A response factor less than 
zero will always introduce damping into a system. 



Figure 4. - Frequency response at three radial 
positions. Drop conditions: radius, 

50 microns; velocity, 100 feet per second; 
temperature, 650° R. Gas conditions*, pres- 
sure, 300 pounds per square inch; final 
velocity, 800 feet per second; oscillation 
peak-to-peak amplitude, 0. 2 of average gas 
pressure. 
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The frequency for a zero response factor, therefore, divides the frequency response into 
regions of driving and damping and may be considered a critical frequency. It is shown 
to be about 2500 cps for the average response curve at the reference boundary conditions. 
This frequency provides a point of comparison in evaluating the effect of boundary condi- 
tions on frequency response. 


Drop Size 

The effect of increasing drop radius from 50 to 500 microns is shown in figure 5. 
The critical frequency, as previously defined, decreases with an increase in drop size. 
The response factor approaches the same constant negative value independent of size. 
The effect of drop size is principally a frequency displacement of the response curve. A 
peak response at about 100 cps is shown for 150-micron drops, and presumably, a peak 
for 500-micron drops exists at a lower frequency. 

Combustor Pressure 


A change in mean combustor pressure affects the response, as shown in figure 6. 
Increasing pressure increases the critical frequency and also gives a higher peak re- 
sponse. This reflects the increase in both acceleration forces and thermal driving force 



Figure 5. - Effect of initial drop radius. Drop 
conditions: velocity, 100 feet per second; 
temperature, 650° R. Gas conditions: pres- 
sure, 300 pounds per square inch; final 
velocity, 800 feet per second; oscillation 
peak-to-peak amplitude, 0.2 of average gas 
pressure. 



Figure 6. - Effect of combustor pressure. 

Drop conditions: radius, 50 microns; veloc- 
ity, 100 feet per second; temperature, 650° R 
Gas conditions: final velocity, 800 feet per 
second; oscillation peak-to-peak amplitude, 

0. 2 of average gas pressure. 
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with pressure and also the change in properties as the drop heats to higher temperatures 
at higher pressures. It also may reflect a limitation of the analytical model. When drop 
temperatures approach the thermodynamic critical temperature, vaporization rate is ac- 
celerated, but the accuracy of the analytical expressions becomes uncertain (ref. 6). Os- 
cillations in pressure, superimposed on the high chamber pressure, may have caused the 
drops to approach this region of accelerated vaporization and, thus, given an uncertain 
peak value. 


Final Gas Velocity 

Response curves for final gas velocities of 400, 800, and 2400 feet per second are 
shown in figure 7. Final gas velocities of this magnitude are reached in combustors with 
contraction ratios of about 6, 3, and 1. 5, respectively. Increasing final gas velocity in- 
creases the critical frequency, but it has little effect on the shape of the response curve. 


Initial Drop Velocity 


Initial drop velocity has a negligible effect on frequency response, as shown in fig- 
ure 8 for injection velocities of 50, 100, and 200 feet per second. An increase in initial 
drop velocity decreases the distance for complete vaporization in steady combustion 
(ref. 3), but it appears to have little effect on the time- dependent properties evaluated in 
this study. 



Figure 7. - Effect of final gas velocity. Drop 
conditions: radius, 50 microns; velocity, 
100 feet per second; temperature, 650° R. 
Gas conditions: pressure, 300 pounds per 
square inch; oscillation peak-to-peak 
amplitude, 0.2 of average gas pressure. 



Figure 8. - Effect of initial drop velocity. 
Drop conditions: radius, 50 microns; 
temperature, 650° R. Gas conditions: 
pressure, 300 pounds per square inch; 
final velocity, 800 feet per second; 
oscillation peaK-to-peak amplitude, 0.2 
of average gas pressure. 
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Figure 9. - Effect of initial drop temperature 
Drop conditions: radius, 50 microns; 
velocity, 100 feet per second. Gas condi- 
tions: pressure, 300 pounds per square 
inch; final velocity, 800 feet per second; 
oscillation peak-to-peak amplitude, 0.2 of 
average gas pressure. 



Figure 10. - Effect of peak-to-peak pressure 
amplitude ratio. Drop conditions: 
radius, 50 microns; velocity, 100 feet per 
second; temperature, 650° R. Gas condi- 
tions: pressure, 300 pounds per square 
inch; final velocity, 800 feet per second. 


Initial Drop Temperature 

The effect of changes in initial drop temperature from 500° to 800° R is shown in 
figure 9. Drop injection temperature affects the magnitude of the response factor only at 
low frequencies, but it has no appreciable effect on the critical frequency. Peak positive 
response factors increase with a decrease in initial temperature. 


Pressure Amplitude 

Figure 10 shows that peak-to-peak pressure amplitude ratios of 0. 1, 0. 2, and 0. 4 
have a small effect on the magnitude of the response factor; the critical frequency, how- 
ever, increases with an increase in wave amplitude. This characteristic indicates a 
changing frequency response during the buildup period of pressure oscillations in an un- 
stable combustor. 


Propellant Properties 

The frequency response for liquid oxygen drops vaporizing in combustion gases com- 
posed of stoichiometric reaction products with hydrogen was evaluated at the reference 
boundary conditions used for heptane vaporization but with an injection temperature of 
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Response factor, N, dimensionless 



Figure 11. - Effect of propellant properties. Drop 
conditions: radius, 50 microns; velocity, 100 feet per 
second; heptane temperature, 650° R; oxygen temper- 
ature, 140° R; radial position, 0.707 of combustor 
radius. Gas conditions: pressure, 300 pounds per 
square inch; final velocity, 800 feet per second; oscil- 
lation peak-to-peak amplitude, 0. 2 of average gas pres- 
sure. 


transformation in frequency, as shown in fi 
giving this correlation of results is 


140° R. This exploration of the effect of pro- 
pellant properties is shown in figure 11. The 
significant differences in the frequency re- 
sponse for oxygen compared with heptane are 
the larger peak values of the response factor 
and the absence of a region of negative re- 
sponse factor. This comparison shows that 
propellant properties can have a significant 
effect on frequency response of a vaporiza- 
tion process. 

Correlation of Results 

The frequency response curves for hep- 
tane vaporization were combined to obtain a 
single curve representing the results for the 
range of boundary conditions evaluated in 
this study. The curves were combined by a 
;ure 12 (p. 18). The transformed frequency 



At the critical frequency, F is 2 500 cps. The factor is normalized to the reference set 
of boundary conditions to facilitate numerical evaluations. 

The complete range of boundary conditions correlates with the frequency factor to 
within a relatively small deviation. Figure 13 shows, by the shaded area, the range of 
deviation within which all the results fall. The greatest deviation exists at the peak val- 
ues. A single curve, however, is adequate to represent the frequency response to within 
the accuracy expected for the analytical model and method of evaluation. 

The curve is applicable only to heptane vaporization with traveling transverse acous- 
tic oscillations. Other modes, principally standing, prescribe relations between pres- 
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Response factor, N, dimensionless Response factor, N, dimensionless 



Transformed frequency, F, cps 


Figure 12. - Correlation of representative heptane results. Initial axial drop velocity, 
100 feet per second; initial drop temperature, 650° R. 



Transformed frequency, F, cps 


Figure 13. - Deviation from single curve in correlation of all results. 
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sure and velocity oscillations which may differ substantially. An effect on frequency re- 
sponse is expected. 


ANALYSIS OF RESULTS 

The effect of boundary conditions on the frequency response of the vaporization pro- 
cess is, in most instances, characterized by a frequency displacement of the response 
curve. This forms the basis of the correlation in figure 12. An alternative treatment of 
the data is to relate the response factor to the characteristic times of the processes in- 
volved. These times are the lifetimes (or half- lifetimes) of the vaporizing drop and the 
wavetime of the oscillation. Figure 14 shows the response factor as a function of the 
ratio of drop half-lifetime to wave or cycle time. The correlation obtained is qualitative- 
ly equivalent to that of figure 12. Peak and zero values of the response factor occur at 



Drop half-lifetime/Cycle time, dimensionless 


Figure 14. - Correlation of representative results with dimensionless time. 
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Figure 15. - Effect of restrictions on drop be- 
havior. Drop conditions; radius, 

50 microns; velocity, 100 feet per second; 
temperature, 650° R for reference case, 

840° R for all others; radial position, 0. 707 
of combustor radius. Gas conditions: pres- 
sure, 300 pounds per square inch; final 
velocity, 800 feet per second; oscillation 
peak-to-peak amplitude, 0.2 of average gas 
pressure. 


dimensionless times of about 0. 4 and 1. 4, respec- 
tively. This means that the response factor is 
near a maximum when the entire drop lifetime 
equals the oscillation period and negative when the 
lifetime is greater than three times the period. 

In an attempt to relate specific characteristics 
of the response curve to their causes, detailed time 
variations in vaporization parameters were exam- 
ined in several drop histories. In addition, calcu- 
lations were made in which the behavior of the drop 
was restricted by holding terms in the vaporization 
equations constant. The results are shown in fig- 
ure 15, with an unrestricted reference case for 
comparison. 

The response factor (fig. 14) is negative when 
vaporization time is large compared with cycle 
time, and thus, many oscillations occur in the 
vaporization history. In this region, the thermal 
inertia of the drop is high relative to the oscilla- 
tion frequency, and the mean drop temperature 
rises from its injection value to the equilibrium 
temperature with only a small oscillatory compo- 
nent. The mass transfer can be related to the 
pressure parameter 


w ~ 



(17) 


where the vapor pressure at the drop surface P L is constant for a constant temperature. 
Under these conditions, this analytic form predicts a decrease in vaporization rate with 
an increase in pressure, and it is the cause of negative response when drop lifetime is 
large. This is illustrated in figure 15, curve D, for computations in which drop tempera- 
ture is held constant at the equilibrium value. The response factor is negative for all 
frequencies. 

As drop lifetime approaches the period of oscillation, the drop temperature can re- 
spond to the acoustic oscillations. This favors mass transfer in phase with pressure and 
results in a peak positive response factor. 

The elimination of the heating period by the injection of the drop at the equilibrium 
temperature (fig. 15, curve C) is amenable to the interpretation just given. At high fre- 
quencies, the heating period is relatively unimportant, because the drop temperature 
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changes very little during one cycle of oscillation throughout its lifetime. By the previ- 
ous consideration,, the response factor is negative for this condition and nearly identical 
to a drop with a heating period (curve A). In contrast, the heating period affects the re- 
sponse of a drop with a lifetime nearly equal to the cycle time. Gains during this heating 
period are larger than at equilibrium temperature and contribute to an increase in over- 
all gain. The peak positive response factor, therefore, is largest for a drop with a heat- 
ing period. The effect of initial temperature on response factor (fig. 8) results from a 
change in the heating period. 

When drop lifetime is very small compared with cycle time, the response factor 
tends toward zero. In this region, the pressure and transverse gas velocities are essen- 
tially constant during the drop lifetime. In the limit of very small drops, the propellant 
is consumed as fast as it is injected regardless of chamber conditions. This gives a re- 
sponse factor of zero. 

Variations of p AV in the Reynolds number also cause mass transfer perturbations. 
Velocity difference appears as a rectified function with almost identical curves over both 
half cycles of oscillation. Density, however, is in phase with pressure, and the product 
p AV produces only in-phase perturbations between mass transfer and pressure. This 
gives a positive response factor, as shown by curve B, where only the Reynolds number 
in the vaporization equation was allowed to vary. 

The drop acquires transverse velocity components because of acceleration by trans- 
verse gas velocities. The phase between the transverse gas and drop velocities affects 
the phase and magnitude of AV, especially for small, easily accelerated drops. To test 
the significance of this effect, curves B, C, and D of figure .15 were recalculated with 
transverse drop velocity held at zero. No appreciable change occurred in the shapes of 
the curves for the conditions examined. 

APPLICATION OF RESULTS 

The application of the response characteristics of the vaporization process to com- 
bustor design may be approached in several ways. According to the correlation parame- 
ter of figure 12, a stable design with the response factor less than zero requires that 
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The frequency of oscillation is related to the combustor radius and the speed of sound by 
the relation (ref. 4) 


f = 


0.293 — 


*w 


(19) 


If an average speed of sound of 42 50 feet per second is assumed for hydrocarbon and oxy- 
gen combustion products, the criterion for a zero response factor becomes 



As a scaling parameter for R w and R d Q , equation (20) states that the drop-size ratio 
must increase by the 2/3 power of the combustor radius to hold the response factor con- 
stant. When combustor thrust level increases by the square of the combustor radius, 
the drop size must increase by the 1/3 power of the thrust-level ratio. The other param- 
eters in equation (20) may be examined for scaling in a similar manner. 

In some applications, the combustor length required to vaporize or burn the propel- 
lants may be a design parameter that is more convenient to use than the vaporization pa- 
ramater in equation (20). An approximate stability criterion involving the vaporization 
length can be formulated from the dimensionless time correlation (fig. 14, p. 19). This 
correlation predicts a response factor less than zero when 

tso 

— >1.4 (21) 

t 

The dimensionless time can be transformed into a parameter ratio involving vaporization 
length from the following manner. 

The oscillation period is related to the combustor radius, as in equation (19), 

R 

t^ = 3. 41 — (22) 

~ a 

The length to vaporize 50 percent of the drop weight L^q is equal to the product of 
the drop half- lifetime and an average drop velocity, so that 
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(23) 


t 



The drop velocity (from a generalization of vaporization histories) is assumed to in- 
crease linearly with combustor length to a final value equal to 80 percent of the final gas 
velocity. This gives an average value of 20 percent of the final gas velocity during the 
half- lifetime of the drop. Final gas velocity or Mach number is related to the combustor 
contraction ratio si- This relation leads to the approximate equation 


t 


50 


7.7 




(24) 


When equations (22) and (24) are combined, the expression for dimensionless time 
is 


50 


= 2.26 



(25) 


When combined with equation (21), this criteria for a negative response factor is 


/VU> 0 .6 

V V 


(26) 


Equation (26) states that Lgg/R^ must remain constant when scaling at constant 
contraction ratios and increase with any decrease in contraction ratio. This conflicts 
with general practice. Large combustors of high thrust are usually designed with 
smaller contraction ratios and aspect (length-to-diameter) ratios than small combustors 
of low thrust. This prescribes a lower value of L^q/R^^ for high-thrust combustors 
than the low-thrust combustors when the performance-length relation is maintained simi- 
lar. In practice, therefore, the high- thrust combustors may be expected to have smaller 
values of (L^q/R w )j/ than low-thrust combustors. This could lead to positive response 
factors and greater instances of instability in high-thrust combustors when the value of 
(L^q/R w )^ is near the critical value of 0. 6. 

The value of for some hydrocarbon and oxygen combustors of high 

thrust, which experienced transverse instability, can be computed from the data in refer- 
ence 8. These are shown in figure 16. Stable and unstable configurations are plotted on 
the response curve, to show both (L^q/R^.^ and the expected value of the response 
factor. Figure 16 shows that these high-thrust combustors are in the region where in- 
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Illl 1 1 1 I I till 1 1 I I I I 
10 1 .1 
Burning rate parameter, dimensionless 

Figure 16. - Response factor related to combustor variables. 

stability is most probable for a vaporization- limited combustion process. 

The correspondence between (L^/R^) d and the burning rate parameter used 
in the stability analyses of references 1, 2, and 8, is also shown in figure 16. They are 
related by the expression 


<27) 

y mR w Vv 

when m, the weight fraction vaporized per unit combustor length, is equated to 
0. 5/L5Q. Relating the response factor to Sf, as in figure 16, adds a new dimension to 
interpreting the stability criteria presented in references 1 and 2. The response charac- 
teristics of the vaporization model used in these studies depended only on Reynolds num- 
ber variations (fig. 15, curve D) and did not give a negative response region. 


COMPARISON WITH CHARACTERISTIC TIME AND INTERACTION INDEX 

It is postulated in reference 7 that the combustion process in a rocket combustor 
may be characterized by two parameters, a characteristic time r and an interaction in- 
dex n. System stability is analyzed by using these parameters, and stability limits are 
analytically predicted for a variety of system boundary conditions. Such a characteriza- 
tion of the combustion process can be used to approximate the response of the vaporiza- 
tion process over a limited range of frequencies. 
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The expression for the fractional change in burning rate in terms of the characteris- 
tic time and interaction index is given in reference 6 as 


w 



w 




(t - t) 


where burning rate and vaporization rate are considered to be identical. If a sinusoidal 
pressure oscillation with time is assumed, the response factor (eq. (13)) for this func- 
tional form can be analytically evaluated and is given by 

N = n(l - cos 27rfT) (27) 

This response function has a peak value of 2n at a frequency fp ea ^ equal to 1/2 t 
cps. The peak amplitude of the response factor obtained from the average curve through 
the correlated vaporization results in figure 13 (p. 18) was about 0.74, which corre- 
sponds to an n equal to 0. 37. This peak occurred at a transformed frequency value 
Fpeak of 750 C PS. This value may be related to f peak by the frequency factor (eq. (16)) 
to give the following expression for characteristic time: 




Figure 17. - Comparison of frequency response of vaporization 
process with function containing characteristic time and interaction 
index. 


Setting n equal to 0.37, using equa- 
tion (28) for t, and relating f to F 
again by equation (16), give the re- 
sponse curve shown in figure 17. It is 
seen that this response function ap- 
proximates the low frequency portion 
of the response curve of the vaporiza- 
tion process. 

Another general expression for 
characteristic time may be obtained 
from the correlation parameter dimen- 
sionless time t^/t^. Figure 14 
(p. 19) shows that the response factor 
peaks at a dimensionless time equal to 
0. 4. By the inverse relation between 
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wave time and frequency, 


T 



(29) 


An interaction index of 0. 37 is in general agreement with what has been postulated 
for the combustion process in liquid propellant combustors. A decrease in characteris- 
tic time with a decrease in drop size and an increase in combustor pressure as given by 
equation (28) or, more generally, a characteristic time proportional to the time to burn 
has been anticipated in liquid propellant combustors. 

It may be concluded from this discussion that the response of the vaporization pro- 
cess may be partially approximated by a characteristic time and an interaction index and 
used to establish stability limits of systems previously analyzed with these parameters. 


SUMMARY OF RESULTS 

The frequency response of the vaporization process was evaluated from an analysis 
of drop vaporization in a combustor with superimposed traveling transverse acoustic os- 
cillations. The response is expressed as the ratio of the percentage oscillation in vapor- 
ization rate in phase with pressure to the percentage oscillation in pressure. The follow- 
ing results were obtained: 

1. A typical response curve attained a peak value at a specific frequency. It de- 
creased to zero at lower frequencies and to a constant negative value at higher frequen- 
cies. 

2. Frequency response characteristics were related to a dimensionless time, which 
is the ratio of the time required to vaporize 50 percent of the drop mass to the period of 
the oscillation. Peak response occurred at a dimensionless time of 0. 4. Negative gain 
occurs at dimensionless times greater than 1. 4. 

3. Frequency response curves for heptane were correlated by a transformed fre- 
quency factor. The factor at the critical condition for zero response factor is 



where q is the drop radius, P c is the combustor pressure, Up is the final axial 


26 



gas velocity, and AP/P c is the pressure amplitude ratio. Two other parameters exam- 
ined, initial drop velocity and temperature, had a small effect on the critical frequency 
dividing regions of positive and negative gain. 

4. The characteristics of the frequency response curves were physically related to 
the temperature of the drop. A large oscillating component in drop temperature gave a 
positive response factor, while a drop without temperature oscillations gave a negative 
response factor. Reynolds number variations always contributed to a positive response 
factor. 

5. Calculations for oxygen vaporization showed propellant properties to affect fre- 
quency response. The peak value and corresponding frequency required for the peak 
value of the response factor were significantly higher than for heptane. No region of 
negative gain occurred for oxygen. 

6. A comparison of results with combustor design practices implied that instability 
is more probable in high-thrust combustors than in low- thrust combustors. 

7. A comparison of analytical and experimental results showed that unstable combus- 
tor designs are in the region where instability is most probable for a vaporization- limited 
combustion process. 

8. The response of the vaporization process may be partially approximated by a 
characteristic time and interaction index previously used to characterize the combustion 
process in system stability analysis. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 14, 1966. 
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